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Charge-Density- Wave Formation in the Doped Two-Leg Extended Hubbard Ladder 
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We investigate electronic properties of the doped two-leg Hubbard ladder with both the onsite and the 
nearest-neighbor Coulomb repulsions, by using the the weak-coupling renormalization-group method. It 
is shown that, for strong nearest-neighbor repulsions, the charge-density-wave state coexisting with the 
p-density-wave state becomes dominant fluctuation where spins form intrachain singlets. By increasing 
doping rate, we have also shown that the effects of the nearest-neighbor repulsions are reduced and the 
system exhibits a quantum phase transition into the J-wave-like (or rung-singlet) superconducting state. 
We derive the effective fermion theory which describes the critical properties of the transition point with 
the gapless excitation of magnon. The phase diagram of the two-leg ladder compound, Sri4_j:Caj;Cu2404i , 
is discussed. 

KEYWORDS: doped Hubbard ladder, intersite Coulomb repulsion, spin gap, charge density wave, su- 
perconductivity 



Electronic properties on ladder systems have been studied 
' intensively both theoretically and experimentally, since the 
superconducting (SC) state was discovered in the self-doped 
two-leg ladder material Sri4_^Ca;tCu2404i with > 12 under 
. pressure over 3 GPa.^'^ The substitution of Ca for Sr changes 
' effectively the hole-doping rate in the ladder Cu sites, where 
. the rate varies monotonically from 0.07 to 0.25 with increas- 
' ing X from to 12.^^ A characteristic feature is the presence 
, of a gap in magnetic excitations at temperature much higher 
' than the SC transition temperature.'^'^ Besides the SC state, 
recent experimental studies have focused on the charge dy- 
' namics in the slightly doped materials and verified collective 
, modes from the sliding of the charge-density-wave (CDW) 

■ developed on ladder sites.^"^ A global phase diagram is ob- 
, tained on the plane of x and temperature showing that the hole 
' doping suppresses the CDW state followed by the insulating 
, state without the CDW order, and then the high doping leads 

■ to the SC state under pressure.^ Quite recently the CDW col- 
lective modes are also suggested in the highly doped material 
Sr2Cai2Cu2404i.^'' Therefore it is of particular interest to in- 
vestigate the competition between the SC state and the CDW 
state in doped ladder systems. 

\ From a theoretical point of view, the origin of the spin gap 
in the ladder compounds seems to be explained successfully 
for both the undoped^'' and doped'^ cases, and it is known 

■ that the c/-wave-like SC (SCd) state appears in doped ladder 
systems.'^ In addition, the charge-ordered state is also sug- 
gested when intersite interactions are included.'^ Further the 
competition between the SCd state and the charge-ordered or 
charge-density-wave (CDW) state has been examined.'^' 
However, critical behavior is not yet fully understood. In the 
present paper, the possible scenario of the instability of the 
CDW state and the competition between the CDW state and 
the SC state are proposed in the doped two-leg ladder of the 
extended Hubbard model (EHM) with nearest-neighbor re- 
pulsive interactions by extending the previous analytical cal- 
culations'^' to analyze the critical behavior in more detail. 

We consider the two-leg EHM given by // = Ho + Hmt- The 
first term describes the hopping energies along and between 



legs: 

^0 = -t\\ L(4/.aC;+l,/,a + H.C.) 

-'^L^'^lLa'^JXa + a-^-), (1) 

where Cjj a annihilates an electron of spin (7(=t, j) on rung 
j and leg l{— 1,2). The Hamiltonian Hint denotes interactions 
between electrons: 

^^int = 1^ n;,/,T n;,/,i + ^11 1^ n j./ n j+ 1 ,/ + V± XI ";M ' (2) 

jJ j-I j 

where U represents on-site repulsion and V\\ (V±) repre- 
sents intrachain (interchain) nearest-neighbor repulsion with 
njj,a (^],i.a^ i,i.c ^"'l = "j,'.T The Hq term is 

diagonalized by using the Fourier transform of Ca{k) where 
k — ,fcx) with k±=0 or n. The energy dispersion is given 
by e(k) = — 2f|| cos^n — t±cosk±. Here we consider the case 
with finite hole doping 8 satisfying tj_ < 2f|| cos^ ^5, in which 
both the bonding {k± — 0) and the antibonding {k± ~ n) en- 
ergy bands are partially filled and the Fermi points are lo- 
cated at kf.Q ~ f (1 — 5) +A and at kp^Tc = f (1 — 5) — A with 
X = sin^' [fx/ (2f|| cos -15)] . We examine the case of small 
5 by neglecting the differences in the Fermi velocities of the 
bonding/antibonding band, i.e., Vf o = Vf^-,i{= vp). 

Following the standard weak-coupling approach (g-ology), 
the Unearized kinetic energy is given by Hq = Y.k.p.a {p^\\ ~ 
kfii^)c^p(j{k)cp„{k), where the index p — +/— denotes 
the right-Zleft-moving electron. By introducing field oper- 
ators by Wp^a.dx) - i"'/^n||e'*^«%,a(^||,^±) with C = 
+ (— ) for k± — 0{7t) and L being the system size, the in- 
teractions near the Fermi points are rewritten as Hint = 
(1/4) /<ixLp,aLf,=±'^nt, where is given by 

^1^2)11 "^/L.Ci ^+{-)P,<y,QA "^-(+)p,<T.f3 

and CT =t (i) for a =[ {]), e = Ci Cs and e = Ci C2- The sum- 
mation of the band index i^, {i — 1, ... ,4) is taken under the 
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condition Ci C2 ^3 Ct = + 1 • Th^ coupling constants and gf^^ 
with / = 1 (2) corresponding to the backward (forward) scat- 
tering are given by = {leV± +mi^eV\\) and gfl = {U + 
/gVx +'Mi\£^|) where /± = ±1, mi_+ = — 2cos7r5 cos2A, 
mi _ = — 2cos7r5, OT2.+ = +2, and mj,- = +2cos2A. We ne- 
glect the umklapp scattering processes which become irrele- 
vant for finite doping case and also neglect forward scattering 
processes which do not yield qualitative changes in the sys- 



tem 



16 



As possible states, we consider the singlet <i-wave super- 
conducting (SCd) state, the CDW state, and the /?-density- 
wave (PDW) state. The PDW state corresponds to the spin- 
Peierls state in the limit of 5 ^ 0. The order parameter 
of the SCd state is given by Oscd = -^"^ Lj(c/,i,t c,-,2,i — 
C;,i,j,cy^2,T). while those of the density waves are Oa = 
N-''Lk,afA{k)cUk)c^{k + Q), with Q= (7r(l-5),7r), 
/cDW = 1 and /pDW = sinA:||. These operators are rewrit- 
ten in terms of bosonic phase fields by applying the 
AbeUan bosonization method.^^'^^ The field operators of the 
right- and left-moving electrons are written as Wp,a,i;{^) = 
ri^^{27ta)^^^^ exp[ipkf i^^x + ip(Pp^^^{x)] where s = + for 
a =t and s = — for a =J, and the fields satisfy the 
commutation relations: [(Pp^s,^{^)t9p.s' .^'{^)] = ip7tsgn{x — 
x') 5, 5/ 5^^/ and [9+..s.f . <P-..v'.f'] = iT^^s/ The Klein fac- 
tors r](j are introduced in order to retain the correct anticom- 
mutation relations.^'* For calculating physical quantities, the 
field (Pp^s i^ is replaced by new bosonic fields: 0vr = (0vr + 

and 0vr = {^tr - ^vr) whcre = ((^>^+ + C(^>^_ + 

^i'a+ + ^C^a-) with p = ±, s = ±, and ^ = ±. The phase 
fields (j)p± and (pa± represent charge and spin fluctuations, re- 
spectively and the suffices ± refers to the even and odd sec- 
tors. They satisfy [^vr('t), 0vv(^)] = —i7t&{—x+x')dr/ with 
&{x) being the Heaviside step function. In terms of <j)vr and 
Ovr, the order parameters O = f dxff are given by 



^S,Cd °= e P+ COS 0p_ COS COS 

— je'^P^ sin Op- sin sin 0cj_ , 
^CDW °= COS0P+ sin0p_ cos(^cy+ cosQ^-, 

— sin^p-i- cos0p_ sin^CT+ sinQ^- 
^PDW °^ cos(^p+ cos0p_ sin(^cy+ sinQ^-) 

-h sin (^p_|_ sin 0p_ cos <pc+ cos Qa- 



(4a) 



(4b) 



(4c) 



We can also rewrite the Hamiltonian in terms of the bosonic 
phase variables. In Eq. (3), the phase field (^p„ appears in 
the form cos(2(^p_ +AXx). Since we can safely assume that 
t±_ is a relevant perturbation for t±_ being not so small, ^^''^ 
the term with cos(20p_ + AXx) would become irrelevant, 
and thus we discard it in the following. We also neglect the 
cos20CT- cos20CT- term since its scaUng dimension is larger 
than 2. Then our Hamiltonian reduces toH = f dxJif with 



E 

r=± 



■P=± 



8pr 



VF 

1 



+ ^ E [ E Wrf - 1^ [dA^r) {dAar) 

Jl r 1- >-„ 1_ JI,VF 



2;rV 



+ g-^;^ cos20CT_) cos20p_ 

{gs+,s- COS20O-+ COS20(J- 

+ g^^ — cos2(pa+ COs20fj-), 



(5) 



where the coupling constants of the harmonic terms are given 
by gpia)r = i:s=±/p%y(4r + Hstl - gll)/2 with r = ±, 



= 1, = e, /J+ = -1 and /J_ = -e. The coupUng 

constants of the nonlinear terms are g— = —^1^, g— v- = 

-82±' 8—,— = {82^ ~ 8i\\^)' 8s+,s- = 8i±' 8s+,~ = 8i± ■ 
These nine coupUng constants are not independent since 
the global spin-rotation SU(2) symmetry leads to'"* {8a+ + 

8a- - 8s+,s-) = {8a+ - 8a- - 8s+^) = (^^,*+ - 8-^,s- - 
8-^j=) = 0. We choose following six coupling constants: 



8-d=,st = +U-V±- 2V\\ cos 7t5, , 
8—,ss= +f/-Vj_-h2V^|(cos;r5-F2cos2A), 
gp+ = +U + 2V± + Vii[4 + cos;r5(l -I- cos2A )] , 
gp- = —Vj_ — V^|Cos;r5(l — cos2A), 
= — ^|Cos;r5(l -|-cos2A), 
gc- = -l-yi-l-y||Cos7r5(l — cos2A). 



where g-^ 



■,s+ 



and g-^ 



(- 



+ 1 



(6a) 
(6b) 
(6c) 
(6d) 
(6e) 
(6f) 



The present model and the above treatment are quite sim- 
ilar to those in Ref. 15. However, the application of the 
renormalization-group (RG) method to Eq. (5) is compUcated 
to estimate excitation gaps of spin modes properly. Therefore, 
we fermionize the spin part of Eq. (5)'"* by introducing spin- 
less fermion fields y/±A^) = Tjr(27ra)~^/^exp[±j20^r(x)] 
where r = ± and {77^, rjr'} = ^^r.r'- By using the SU(2) con- 
straints and the Majorana fermions ^" (n = 1 ~ 4), the equa- 
tion (5) is rewritten as 

^ = 'iL\Ud^Prf + ^{d.<^;r)id.^pr) 
Jl r p nvp 



.VF 



4 c4 



~ 2^ (^-.^^^+ • +S-,ss^t ■ cos20p_,(7) 

where VAp,+ = {^l + i^j) / ^/2, Xj/p,. = {^^^ + i^'^) / V2, and 
= {^p,^p,^p). Thus the effective theory for the spin sector 
becomes 0(3) XZ2 symmetric, as seen in the isotropic Heisen- 
berg" and half-filled Hubbard ladder. 

We investigate the low-energy behavior by using the per- 
turbative RG method with the lattice constant a ae^'. Fol- 
lowing six scaUng equations are obtained: 



dl 
d_ 
Jl 



3 2 1 
:Go- = —-G—,, -G 



4 ~ c—,ss^ 
a+ = -G(j+ - Gfy_ - -U- 



d 



-2Ga+ Ga- 



_rz — n — 

2 c— ,sf ^c—ss' 



(8a) 
(8b) 
(8c) 



dl 



G^=^st = -Gp- G^gt ~ '^Ga+G^^ - Ga- (8d) 
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—G — 



(8e) 



and dGp+/dl = where G(0) = g/ilnvp). We noted 
that these RG equations can be also derived directly 
from Eq. (5). We analyze the RG equations numeri- 
cally for U > Q, V\\ > Q and Vl > 0. For small V^^/U 
and V\\/U, the limiling behavior of RG equations is 



given by {G*p_,G%+,Gl_,G*-^^^,G*-^^ 
which corresponds to {gi- 



',s+ ' ^c- 



* 

■:Ss+,S- 



,8 



s+,s- 



:) = 



(— , — ,0,— ,0) in Eq. (5). The relevant behavior of cou- 
pling constants implies that the phases are locked in or- 
der to minimize the cosine potential in Eq. (5). The posi- 
tions of phase locking and the corresponding ground states 
are summarized in Table 1. Since the 6a- field is conju- 
gate to these two fields cannot be locked at the same 
time. From Eq. (4), the nonvanishing order parameter is 
^scd- Since the correlation function of the operator e'^P+ 
exhibits power-law behavior, we obtain that the SCd fluc- 
tuation becomes quasi-long-range ordered (quasi-LRO) in 
this case. We note that the SCii state moves to the D-Mott 
or D'-Mott state in the hmit of 5 ^ 0.'"* For large V±/U 
and Vii/f/, the limiting behavior of RG equations is now 
given by {G;_,GU,G*a-,G*-^,,G*-J = (--,+,-,+), 
corresponding to ^._,^*^_) = 

0, -|-, 0, — ). In this case, the dominant order parameters are 
given by ^cdw and ^pdw both of which lead to the quasi- 
LRO with the same exponent of the correlation functions. We 
call this coexisting state the CDWh-PDW state. 

In order to analyze the properties near the critical point 
of the transition between the SCd state and the CDWh-PDW 
state, we restrict ourselves to the case where the mass of the 
charge mode (p— ) is larger than those of the spin modes 
(a±). The Op- field is locked by the cosine potential be- 
low the mass scale of the charge mode mp_. By replacing 
cos20p_ with its average value = (cos20p_) in Eq. (5), 
the effective low-energy Hamiltonian for the spin degrees of 
freedom is obtained as''*- 



2^ 



.VF 



2 



(9) 



where Wj and represent bare masses of the Majo- 
rana triplet and singlet sector: = {cp^/2na){U — Vj_ — 
2V^|C0S7r5) and m° = (cp^/27ra)[?7 - Vx + 21^1 (cos 7r5 + 
2cos2A)]. The quantity m'^ (m^) has physical meanings of the 
gap in the magnon (soliton) excitation in the spin modes of 
the ladder.'^ Equation (9) is further analyzed in terms of the 
following scaling equations for coupling constants: 

dGt 
dl 
dGs 
'df 
dGa+ 



= Gt — 2GtGa+ — GsGc- 



—rr =Gs — 3GtGa- 



(10a) 



(10b) 



dl 
dGg. 
dl 



G2 /~t2 /~t2 
a+ - Ga- - G, , 



-2Ga+Ga- — GtGs 



(10c) 



(lOd) 



Table I. Possible states and the corresponding pattern of phase locking 

where / is an integer and the symbol * indicates an unlocked bosonic phase 
field. The signs + and — denote those for renormahzed masses m, and m, 
where \vc have assumed / beine even luiniber. i.e.. c— ', 0. 



State 










mi Ms 


SCd 


¥ 






* 


+ + 


CDW + PDW 




1(^+1) 


* 




+ 



where Gt = /v/r, Gs = nPjvp, and Ga± = ^<T±/27rv/r. The 
couplings Gs and Gt are relevant, while Ga± are marginal. 
In Eq. (10), the Gs term as a function of / increases rapidly 
compared with other G's and becomes relevant at / = 4 cor- 
responding to the energy scale of a gap in the Majorana sin- 
glet mode Wj ~ t\\e^^', where we stop the calculation of Eq. 
(10). The mode remained below the energy scale of is the 
Majorana triplet sector. The effective theory for this mode 
is given by M'^^ = -i\vF{i+ ■ d_4+ - ^_ • d^^-) - imf ^+ ■ 
I- - ka+ (C+ • ^-)' where = Vf [G,(/s) - G^- (/..)] and 
glj^ — 27iVfGcr+{ls)- Then we solve the RG equations given 
by dGt/dl = G, - 2G,G„+ and dGc+/dl = -2Gl+ - Gj 
with the initial conditions Gt{ls) = mf/vp and G(j+{ls) = 
g'lj^/2nvF- We easily find that these RG equations have two 
stable fixed points (G,*,Gct_|_) = {+°°,—'>°) and (— 0°,— °°), 
corresponding to the SCd state and the CDW-hPDW state, 
respectively, where the magnitude of the gap in the Majo- 
rana triplet sector can be estimated from rrit « f||e^''sgn(G^) 
where It is determined by |G((Zr)| = 1 (see Table 1). There are 
also two unstable fixed points (G* , G*^^) = (0, 0) and (0, -00), 
corresponding to the second-order and first-order phase tran- 
sitions,^^ while only the former transition is obtained in the 
present numerical calculation. 

From the numerical integration of the RG equations, we ob- 
tain the ground-state phase diagram shown in Fig. 1 . The SCd 
state (the CDWh-PDW state) is obtained for Vn/U + Vj_/U > 
0.4(< 0.4). The SCd state, on the one hand, is stabilized 
by the on-site repulsive interaction, which segregates up-spin 
from down-spin on the same site and leads to the singlet pair- 
ing on a rung. On the other hand, the CDW+PDW state is 
obtained due to the nearest-neighbor repulsive interactions, 
which induce density wave leading to the singlet state on the 
same site or chain. The effect of is shghtly larger than 
that of Vx although both the intersite interactions have es- 
sentially the same effect of inducing the CDW+PDW state. 
In Fig. 2, the change from the CDW+PDW state to flie SCd 
state is shown with increasing the doping 5(> 0.05). The 
novel aspect of the present paper is the competition induced 
by the doping which reduces the effect of only as shown 
in Eq. (6). In the inset, we show the respective masses esti- 
mated from \ina\ w f|| exp(— (a = t, s, p— ) by noting that 
the corresponding coupling constant \Ga\ becomes of the or- 
der of unity at / = /«. Our system exhibits a second-order 
phase transition and the magnon excitation gap vanishes at 
the quantum critical point (QCP). The critical property for the 
Majorana triplet sector, which differs from that the conven- 
tional Tomonaga-Luttinger liquid, is described by the SU(2)2 
Wess-Zumino-Novikov-Witten model with the central charge 
c = 3/2." 

In the present paper, by applying the weak-coupling RG 
method to the EHM on two-leg ladder, we have shown that 
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Fig. 1. The ground-state phase diagram on the plane of V\\/U and Vj_/t/, 
with J//t|| = 2, 5 = 0.1, and = ?||- 




doping § 



Fig. 2. The doping dependence of the magnon spin gap with U /fy = 2, 
V\\IU = Vl/C/ = 0.25, and t±_ = = 1. In the inset, the doping depen- 
dences of nip-, Ms, and are are shown. 

the doping 8 suppresses the CDW-hPDW quasi-LRO state 
and yields the system to the QCP, and that the SCd quasi- 
LRO state is stabihzed at further doping. Here we discuss 
the experimental results of the two-leg ladder compound 
Sri4_;(CajcCu2404i. The phase diagram of Sri4_j:Ca;cCu2404i 
obtained in Ref. 9 resembles our phase diagram of Fig. 2 if the 
magnitude of gap \m, \ is regarded as the transition tempera- 
ture. On closer look, our phase diagram is contrast to the fea- 
tures of Sri4_;cCajcCu2404i that the resistivity above the tran- 
sition temperatures shows an insulating behavior and there is 
no experimental evidence of the QCP between the CDW state 
and the SC state. In order to explain the phase diagram of 



Sri4_jcCa;tCu2404i, the dimensionality effect and/or the dis- 
order effect has been discussed.^ The quantum critical behav- 
ior would be smeared out by these effects, which are not taken 
into account in the present paper. However it will be still inter- 
esting to examine the competing region in the sense that the 
magnon gap would become extremely small and anomalous 
behavior can be expected at temperatures higher than char- 
acteristic energies of the disorder and the dimensionality. We 
note that the origin of the high temperature insulating phase 
is stiU unknown and the analysis is left for a future study. 
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